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ABSTRACT 

Based on our previous work (jVidotto et aL|[2009aD . we investigate here the effects on the wind and 
niagnetospheric structures of weak-lined T Tauri stars due to a misahgnment between the axis of 
rotation of the star and its magnetic dipole moment vector. In such configuration, the system loses 
the axisymmetry presented in the aligned case, requiring a fully three-dimensional approach. We 
perform three-dimensional numerical magnetohydrodynamic simulations of stellar winds and study 
the effects caused by different model parameters, namely the misalignment angle 9t, the stellar period 
of rotation, the plasma-/?, and the heating index 7. Our simulations take into account the interplay 
between the wind and the stellar magnetic field during the time evolution. The system reaches a 
periodic behavior with the same rotational period of the star. We show that the magnetic field 
lines present an oscillatory pattern. Furthermore, we obtain that by increasing dt, the wind velocity 
increases, especially in the case of strong magnetic field and relatively rapid stellar rotation. Our three- 
dimensional, time-dependent wind models allow us to study the interaction of a magnetized wind with 
a magnetized extra-solar planet. Such interaction gives rise to reconnection, generating electrons that 
propagate along the planet's magnetic field lines and produce electron cyclotron radiation at radio 
wavelengths. The power released in the interaction depends on the planet's magnetic field intensity, 
its orbital radius, and on the stellar wind local characteristics. We find that a close-in Jupiter-like 
planet orbiting at 0.05 AU presents a radio power that is ~ 5 orders of magnitude larger than the 
one observed in Jupiter, which suggests that the stellar wind from a young star has the potential 
to generate strong planetary radio emission that could be detected in the near future with LOFAR. 
This radio power varies according to the phase of rotation of the star. For three selected simulations, 
we find a variation of the radio power of a fact or 1.3 to 3.7, depending on Ot- Moreover, we extend 
the investigation done in iVidotto et al.l (|2009al ) and analyze whether winds from misaligned stellar 
magnetospheres could cause a significant effect on planetary migration. Compared to the aligned case, 
we show that the time-scale t^, for an appreciable radial motion of the planet is shorter for larger 
misalignment angles. While for the aligned case t^ ~ 100 Myr, for a stellar magnetosphere tilted by 
9t = 30°, Tyj ranges from ^ 40 to 70 Myr for a planet located at a radius of 0.05 AU. Further reduction 
on Tw might occur for even larger misalignment angles and/or different wind parameters. 

Subject headings: MHD - magnetic fields ~ methods: numerical - stars: pre-main sequence - stars: 
winds, outflows ~ planets and satellites: general 



1. INTRODUCTION 

T Tauri stars are pre-main sequence low-mass stars 
(0.5 < M/Mq < 2), with a range of spectral types 
from F to M, and radius ^3 — 4 Rq. They are usu- 
ally classified in two categories, depending on their evo- 
lutionary stage. In an earlier stage, they are known as 
classical T Tauri stars (CTTS), surrounded by circum- 
stellar disks. In a later stage, with the dissipation of 
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the accretion disk, they are known as weak-lined T Tauri 
stars (WTTSs). Thanks to spectropolarimetric measure- 
ments, the number of young low-mass stars with detected 
magnetic fiel ds has significantly increased in the past 
decade (iJohns-Krull et a l. 1999; Valenti & Johns-Krul 
2OOI I Johns-Krulll 120071: Donati e t al. 20 07. .2008 .201 
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Yang et al.ll2008t IHussain et al.l l2009: Skc llv et al.bTM 
These detections have suggested that T Tauri stars 
present mean surface field strengths of the order of kC. 
Surface magnetic maps, derived from spectropolarimet- 
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ric data, indicate that the surface fields on T Tauri 
stars are more complex than that of a simple dipole 
and are often misaliKned with the rotational axis of the 
star (iDonati et al.llgOOTl lMl. CT TSs such as BP Ta u 
(jDonatiet a l.ll2008D. ¥2129 Qph (JDonatiet all [2001. 
CR C ha and CV Cha (IHussain et al.l l2nnfd[^ V2247 
Oph (jDonati et al.l 120100 present dipolar and octupolar 
components of the surface magnetic field moment that 
are asymmetric with respect to the rotational axis of 
the star. More recently, the first determination of sur- 
face magnetic maps for a W TTS, V410 Tau, has been 
acquired (Sk ellv et al.l l2010i) . showing that, similar to 
the less evolved CTTSs, V410 Tau also presents a non- 
axisymmetric poloidal field. 

Despite the existing knowledge of the surface magnetic 
fields in young stars, the global structure of the stel- 
lar magnetic field is unknown. Magnetic field extrap- 
olations from surface magnetograms using the potential 
field source surface (PFSS) method has been used to help 
us elucidate th e geometry of the large-scale field aroun d 
T Tauri stars ()Jardine et al.l 120081: IcFegorv et al.ll2008D . 
Such extrapolations, however, neglect the interaction of 
the field with the stellar wind and the temporal evolution 
of the system. Full r nagnetohydrodynani i cs (MHD) nu- 
merica l simulations (jRoussev et al.|[2003t iVidotto et al.l 
l2009al) allow us to study the interplay between the stel- 
lar magnetic field and the wind. In this method, the dy- 
namical interaction of the stellar wind and the magnetic 
field lines is a result of the action of magnetic, thermal, 
gravitational, and inertial forces. MHD simulations can 
be, however, computationally expensive and time con- 
suming. A comparison between PFSS and MHD models 
that used observed surface m agnetic map s as a bound- 
ary condition can be found in iRilev et al.l (|2006( ) in the 
context of the solar wind, showing that PFSS models are 
able to reconstruct the large-scale structure of the solar 
corona when time-dependent changes in the photospheric 
fiux can be neglected, although nonpotential efi^ects can 
have a significant effect on the magnetic structure of the 
corona. 

The accurate determination of the wind properties 
and topology of the magnetic field of a star is nec- 
essary to solve a series of open questions. The rota- 
tional evolution of the star, for example, requires the 
knowledge of the topology of the field, as the shape of 
the magnetic field lines may enhance rotational brak- 
ing caused by t he stellar inagnet i zed winds (Kra ft 1962 
IKawalei< [19881 : IHolzwarthI 120051: IBouvieH l2009t) . Fur- 
thermore, studies of the magnetic interaction between 
a CTTS and its disk require the know ledge of the struc- 
ture of the magnet ic field of the star (jCamenzindl 119901 : 
[Kocnigl 1991; Lon g et all 120081 ). Determining realis- 
tic magnetic field topologies and wind dynamics are 
also key to understand interactions between magnetized 
extra-solar planets and the star such as interactions 
that l ead to planeta ry migr ation (iRomanova fc Lovelace! 
120061: iLovelaceet al . 2008; Vidotto et al.ll2009al) . interac- 
tions between the stellar ma gnetic field and the planetary 
magnetosphere ('Gricfi mcier et al.l 120041: lip et all I2004J: 
Prcussc ct al. 2006, 2007t iStevensI 120051 : iCohen et al.l 
2009; .Fares ct al..,2010,) and also with the planetary at- 



extend the study performed in IVidotto et al.l (|2009aD . 
where the stellar rotation and magnetic moment vec- 
tors were assumed to be parallel. We now consider 
cases where these vectors are not aligned. Some nu- 
merical and theoretical models exist considering the 
case of an oblique magnetic geom etry, mainly ap- 
plied to th e stu dy of pulsars (e.g.. iBogovaloyI 119991 : 
lSpitkovskvll2006HKalapotharakos fc Contopoulo£ll2009D . 
with a few applications to other astrophysical objects 
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1 & WickramasinghG' 1998l : [Romanova et al.ll2003l 
Townscnd ct al. 2007). As a consequence of the 
oblique magnetic geometry, the sy stem loses the axisym- 
metry present in the aligned case ([Vidotto et al.|[2009a) . 
thus requiring a fully three-dimensional (3D) approach. 
We perform here 3D MHD numerical simulations of mag- 
netized stellar winds of WTTSs, by considering at the 
base of the coronal wind a dipolar magnetic field that 
is tilted with respect to the rotational axis of the star. 
Complex, high-order multipoles magnetic field configu- 
ration at the surface may exist, but a dipola r component 
should dominate at larg er distances (e.g., iJohns-KrullI 
120071 : IDonati et al.l [20070 . As the simulation evolves in 
time, the initial field configuration is modified by the 
interaction with the stellar wind, which in turn also is 
modified by the magnetic field geometry. 

The stellar wind of a host star is expected to directly 
influence an orbiting planet and its atmosphere. The 
interaction, for example, of a magnetized wind with a 
magnetized planet can give rise to reconnection of mag- 
netic held lines. Reconnection processes appear in several 
places in the Solar System. E.g., the magnetic field lines 
of the Earth magnetospheric day-side (i.e., the side of the 
Earth that is facing the Sun) are compressed due to the 
interaction with the solar wind, while in the opposite side 
of the Earth's magnetosphere (the night-side), a tear- 
dropp ed-shaped tail is formed (e.g.. lChapman fc Ferrard 
Il930f ). The solar wind interaction with the magnetic 
planets of the Solar System (Earth, Jupiter, Saturn, 
Uranus, and Neptune) accelerates electrons that prop- 
agate along the planets magnetic field lines, prod ucin; 
electr on cyclotron radiation at radio wavelengths (Zar _ 
Il998f l. By analogy to the magnetic planets in the So 
lar System, predictions have been made that extra-solar 
planets shou ld produce cyclotron maser emiss ion (e.g., 
Farreh et all [1999: Zark a et al.l 120011: iLazio e t al. 2001 
FarrcU et al. 2004: Gricfimcier et al. 2004, 2005, 2007a, 13; 
Stevens 2005; Zarka 2007; Jardine & Cameron 2008), if 
they harbor intrinsic magnetic fields, although predic- 
tions also exist for the case of non-magnetized planet s 
([Zarka et al.ll2001l : lZarka|[2007t iGriefimeier eraI[|2007aO . 
Evidence that extrasolar giant planets can be magnetized 
was found by Shkolnik ct al. (2005, 2008), who observed 
modulations of the Ca II H&K lines in phase with plan- 
etary orbital periods on extrasolar planetary systems. 
Such modulations were interpreted as induced activity 
on the stellar chromosphere caused by the interaction 
between the stellar and planetary magnetic fields.^ 

The consideration of a realistic wind is crucial to de- 
termine how the interaction between the stellar wind and 



mosphere ( Griefimeier et al1l2004[ ). 

As a next step towards a more realistic wind and 
magnetic field modeling of WTTSs, in this work we 



^ IPreusse et al.l l|200g ) suggest that a conductor planet moving 
relatively to the stellar wind is also able to generate perturba- 
tions that can t rigger the ch romospheric modulations observed by 
Shkol nik et al.l IJ2005I . I200SI ) without the requirement of a magne- 
tized planet. 
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the magnetosphere of an extrasolar planet occur. Using 
the 3D, time-dependent MHD wind models developed in 
this paper, we investigate the planet-wind interaction. 
Such interaction can give rise to reconnection processes, 
which result in transfer of energy from the stellar wind to 
the planet's magnetosphere. Analogously to the interac- 
tion of Jupiter's magnetosphere with the solar wind, we 
estimate the radio power released from the interaction of 
a close-in giant planet with the wind of its host star. 

We have organized this paper as follows. Sj2] presents 
the MHD numerical model adopted to describe a magne- 
tized stellar wind of WTTSs. In ^ we present the sim- 
ulations made and the results we achieved, along with a 
comparison between wind models with different parame- 
ters. In Sj4l we discuss our stellar wind results, perform- 
ing comparisons with other simpler stellar wind mod- 
els. The interaction between a close-in magnetized giant 
planet and the stellar wind and an estimate of the ra- 
dio power released from this interaction are presented 
in SjSl In !|6l we investigate whether the action of mag- 
netic torques from the stellar wind acting on a close-in 
giant planet is able to cause planetary migration. This 
investig ation extends the one presented in Vidotto et al.l 
(|2009a[ ). where the winds analyzed in that case assumed 
that the axis of rotation of the star and the surface mag- 
netic dipole moment were aligned. In fJ71 we present the 
conclusions. 

2. THE NUMERICAL MODEL 

To perform the simulations, we use the Block Adaptive 
Tree Solar- wind Roe Upwind Scheme (BATS-R-US), a 
3D ideal MHD numerical code developed at the Cen- 
ter for Sp ace Environm ent Modeling at University of 
Michigan (jPowell et al.r i999). BATS-R-US has a block- 
based computational domain, consisting of Cartesian 
blocks of cells that can be adaptively refined for the re- 
gion of intere st. It has been used to simulate the he- 
liosphere (R oussev et al.l 12003] : [Cohcn ct al. 2007 ]), the 
outer - helios phere ()Linde et al.l |l998; Ophcr ct al.i 120031 
2006 , I2007D. coronal ma ss ejections (Manchester et al.l 
2004: iLugaz et all I2005D . the Earth's magnetosphere 
(Ridlev ct al.' '2006") and the magne tosphere of Saturn 
iHanscn ct al. 2005) and Uranus (lT6th et al.l 120041 ) . 
among others. In this work, w e extend the model devel- 
oped in Vidot to~et al.l ()2009a|) to study the wind struc- 
ture of WTTSs, specifically when the magnetic moment 
of the star and the stellar rotational axis are non-parallel. 

BATS-R-US solves the ideal MHD equations, that in 
the conservative form are given by (in cgs units) 
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Fig. 1. — Meridional cut (xz-plane) of the adopted 3D grid in 
the simulations of misaligned magnetospheres, illustrating the re- 
finement at the inner portion of the grid. Immediately around the 
star, the cell resolution is 0.018 ro (ro is the stellar radius), and 
with distance from the star, the grid gets coarser. The coarsest 
resolution is in the outer corners of the grid (not shown above) 
and is 3.75 ro. 

acceleration due to the central body, and e is the total 
energy density given by 



pu 



V 



1 



Btt 



(5) 



where p is the mass density, u the plasma velocity, B the 
magnetic field, p the gas pressure, g the gravitational 



We consider an ideal gas, so p = pkBT /{pmp)^ where feg 
is the Boltzmann constant, T is the temperature, mp is 
the proton mass, p — 0.5 is the mean molecular weight of 
a totally ionized hydrogen gas, and 7 is the ratio of the 
specific heats (or heating parameter). In our simulations, 
we adopt either 7 = 1.1 or 7 = 1.2. 

The adopted grid is Cartesian and the star is placed 
at the origin. The axes x, y, and z extend from —75 rg 
to 75 To, where ro is the stellar radius. For all the cases 
studied, we apply 11 levels of refinement in the simulation 
domain. Figure [T] presents a cut along the meridional 
xz-plane illustrating the refinement at the inner portion 
of the grid. We note that a higher resolution is used 
around the central star. This configuration has a total of 
~ 2.5 X 10^ cells in the simulation domain. The smallest 
cells (closest to the star) have a size of 0.018 ro and the 
maximum cell size is 3.75 ro. 

The star has M* = 0.8 Mq and ro = 2 i?0. The 
grid is initialized with a ID hydrodynamical wind for a 
fully ionized plasma of hydrogen. The solution for Ur (r) 
depends on the choice of the temperature at the base of 
the wind and on 7, and the only physical possible solution 
is the one that bec omes superson ic when passing through 
the critical radius ()Parkedll958| ) . Due to conservation of 
mass of a steady wind (i.e., pUrV^ = constant), we obtain 
the density profile from the radial velocity profile Ur{r). 

The star is considered to be rotating as a solid body 
with a period of rotation Pq = 27r/il, where Q, is the 
angular velocity of the star. Its axis of rotation lies in 
the ^-direction 

n = rjz . (6) 

The surface magnetic moment vector m is tilted with 
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respect to ft at an angle 9t 

m • r2| = mil cos 9 f . 



(7) 



The simulations are initialized with a dipolar magnetic 
field described in spherical coordinates {r, 9, 1^9} by 



Br {to) 



^O'O , 



cos 



' COS 6*4 + sin 6* cos (/J (to) sin ^t) , (8) 



Beita) — — T^ ( - sinflcos^t — - cos 9 cos ip{tQ) sin 9t 



D r-^ \ ^ora sin ^(io) sin 6*4 
B^[to) - —3 



(9) 
(10) 



where Bq is the magnetic field intensity at the magnetic 
poles (where 9 — 9t and r — tq), r is the radial coordi- 
nate, 9 is the co-latitude, and ip is the azimuthal angle 
measured in the equatorial plane. At the initial instant 
io, the vector m is in the xz-plane, tilted by an angle 9t 
in the clockwise direction around the y-axis (Fig. [5^). 

The inner boundary of the system is the base of the 
wind at r = ro, where fixed boundary conditions were 
adopted. The outer boundary has outflow conditions, 
i.e., a zero gradient is set to all the primary variables (u, 
B, p, and p). As the magnetic field is anchored on the 
star, in one stellar rotational period, the surface magnetic 
moment vector m draws a cone, whose central axis is the 
z-axis (Fig. [2b)- Because of that, in the simulations with 
oblique magnetic geometries, the boundary conditions 
are time-dependent and the simulations reach a periodic 
configuration f H3.1l) . 

The MHD solution is evolved in time from the initial 
dipolar configuration for the magnetic field to a fully self- 
consistent non-dipolar solution. The wind interacts with 
the magnetic field lines and deforms the initial dipolar 
configuration of the field. The stellar wind is also modi- 
fied by the magnetic field, i.e., no fixed topologies for the 
magnetic field neither for the wind are assumed. 

3. SIMULATION RESULTS 

Table [T] shows the parameters adopted in the simula- 
tions. Common to all simulations are the magnetic field 
intensity oi Bq — 1 kG at the magnetic poles of the 
star and the temperature Tq = 10^ K at the base of the 
wind. We varied the misalignment angle 9t, the period 
of rotation of the star Pq, the density at the base of the 
wind po (and consequently the plasma-/3 at the base of 
the coronal wind /3o), and the heating index 7. Observa- 
tions of WTTSs show that th ey possess rotation al peri- 
ods ranging from 0.5 to 13 d dMarilh et al."2007^ with a 
distribution peaking at Pq ~ 2 d (Cieza & Balibcr 2007). 
However, we did not consider Pq > 3 d, as the dynamical 
effect of the misalignment is more significant in the cases 
where the star has a low period of rotation. This implies 
that we are in the lo wer range of obs erved rota tional peri- 
ods for WTTSs (e.^.-IBouvier et aLlllQ QS: Ch oi fc HerbsH 
flOOllHerbst et al.l 120021: IRebuU et aril200a ). A descrip- 
tion of the choice of parameters used in t he simulations 
except for 9t) can be found in Section 3 of lVidotto et al.l 
2009af) . The misalignment angle 9t between the rota- 
tional axis of the star and the surface magnetic dipole 
moment vector was chosen to vary from 0° (aligned case) 
to 90°, although surface magnetic maps of T Tauri stars 



(a) 




m(to + ^Po) \ 



(b) 



Fig. 2. — The magnetic field configuration in the simulations of 
misahgned magnetospheres. (a) Orientation of the magnetic field 
fines in the x2-plane at the initial instant to- (b) In oug stellar 
rotational period, the magnetic moment vector m draws a cone, 
whose central axis is the z-axis. Shown above are the magnetic 
moment vector for the initial instant to and half-period later. 6t is 
the angle between the vectors m and Q. 



have shown the existence of smaller angles 



< 



30°, 



IDonati et all 120071 [200l . Simulation TOl considers the 
aligned case. Simulations T02 to T06 represent case with 
different 9t. With respect to our fiducial case T04, we 
varied 7 in case T07, Pq in cases T08 and T09, and (3q 
in case TIO. 

3.1. Time- Dependent Behavior 

Because the star is rotating and the magnetic field is 
asymmetric with respect to the axis of rotation, the sim- 
ulations with an oblique magnetosphere have a periodic 
behavior with the same period of rotation of the star. 
Depending on the physical conditions and grid size, af- 
ter a certain number of stellar rotations, the system has 
relaxed and such periodic behavior is achieved. For ex- 
ample, for a wind expansion velocity of '^ 200 km s~^ 
propagating in a grid size of 75 tq, the time the solution 
will take to relax in the grid is of r^ 75 ro/(200 km s~^) ~ 
6 days. Cases T02, T03, and T04 were run for 10 days. 
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TABLE 1 

The set of simulations. The columns represent, 

respectively: the name of the simulation, the density po 

at the base of the wind, the heating parameter 7, the 

period of rotation of the star pf), the misalignment angle 

6t, AND THE PLASMA-/3 EVALUATED AT THE MAGNETIC POLE /Bq. 

Case TOl represents the aligned case {9t = 0°) and case 

T04 THE fiducial CASE {9t = 30°). 



Name 


PO 


7 


^n 


&t 


1^0 




(g cm-3) 




(d) 


n 




TOl 


1 X IQ-i^ 


1.2 







1/25 


T02 


1 X 10-" 


1.2 




10 


1/25 


T03 


1 X 10-" 


1.2 




20 


1/25 


T04 


1 X 10"" 


1.2 




30 


1/25 


T05 


1 X 10-" 


1.2 




60 


1/25 


T06 


1 X 10-" 


1.2 




90 


1/25 


T07 


1 X 10-" 


1.1 




30 


1/25 


T08 


1 X 10-" 


1.2 


3 


30 


1/25 


T09 


1 X 10-" 


1.2 


0.5 


30 


1/25 


TIO 


2.4 X 10-12 


1.2 


1 


30 


1/100 



T05 for 9 days, T06 for 8 days, T07 for 5 days, T08 and 
T09 for 6 days, and TIO for 7 days. These time intervals 
were sufficient for the solution to relax in the grid. 

Figure |3] illustrates the time-dependent behavior of our 
simulations, where we show meridional cuts of the total 
velocity of the wind for nine instants during one full pe- 
riod of rotation of the star for case T04, our fiducial case. 
The first panel represents a given instant t = ti; the sub- 
sequent panels increase in multiples of 1/8 Pq, until the 
cycle completes at instant t = ti + Pq. The magnetic 
field lines are represented by black lines and the white 
line denotes the contour where the magnetic field changes 
polarity, i.e., when Br ~ 0. It can be seen that the ini- 
tial instant t = ti and the final instant t = ti + Pq oi 
a given stellar rotational period are identical, which ex- 
emplifies the periodic behavior of the simulation. The 
magnetic field has zones of open and closed field lines. 
The zone of closed field lines rotates around the stellar 
equatorial plane (z = 0), as is readily seen if one uses the 
contour-line of i?^ = (white line) as a guide. 

It should be reminded that Br — defines a bi- 
dimensional surface, but because Fig. [3] (as well as sev- 
eral other figures we present later on) is a meridional 
cut, Br = is shown as a contour-line. When there is no 
misalignment {9t = 0), the surface Br = coincides with 
the equatorial plane (z = 0). However, in the case where 
9t ^ 0^ Br = defines a wavy, time-dependent surface 
for an observer that is in an inertial referential frame. 
In the corotating frame of the star, the surface would 
still be wavy, but it will appear as static. In both cases 
{9t = and 9t 7^ 0), such a surface is the locus of points 
at the tip of the closed magnetic field lines. Therefore, 
we expect cusp- like structures (i.e., helmet-streamers) to 
also oscillate with the same rotational period of the star. 
The three-dimensional view of the surface Br = can be 
seen in Fig. S) 

3.2. The effects of the misalignment angle 9t 

By comparing simulations TOl to T04, we can analyze 
the effects caused due to small misalignment between 
the stellar rotation axis and the stellar magnetic dipole 
moment vector in the wind structure and magnetic field 
configuration. 

Figure [5] presents radial velocity color maps for the 



cases T02 {9t = 10°), T03 {9t = 20°), T04 {9t = 30°), 
and the aligned case TOl {9t = 0°) for the entire sim- 
ulation box (shown in the figure is a meridional cut in 
the a;z-plane). The four panels are snapshots taken at 
t ~ 240 h, which is a sufficient time to assure that the pe- 
riodic behavior described in ij3.1l has been achieved. For 
each of the misaligned cases, we note regions of higher ve- 
locities, surrounded by regions of smaller velocities. This 
oscillatory behavior is seen in all the variables of the wind 
and is caused by the precession of the stellar magnetic 
field around the polar axis of the star. We also note that 
the increase in 9t leads to faster winds on average. This 
is caused by the azimuthal derivative (d/dip) terms in 
Eqs. ([U to (HJ, which vanish in the aligned case. In the 
momentum equation [Eq. [2], for instance, these terms 
are: the inertial term 



r smi 



dur 
dip 



due 
dip 



dip 



■ ip 



the pressure gradient in the azimuthal direction 

1 dp „ 
r sin 9 dip 

and the magnetic force (V x B)xB/47r 



(11) 
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•^ 



(13) 

where the first term inside the brackets refers to a mag- 
netic tension and the second term refers to a gradient of 
the magnetic pressure. The terms given by Eq. ([TH) are 
negligible when compared to the total inertial force, as 
well as the terms in Eq. (|12|) when compared to the mag- 
nitude of the total pressure gradient force. The magnetic 
tension and pressure presented in Eq. (|13p contribute 
more significantly to the acceleration of the wind under 
an oblique magnetic field configuration. Furthermore, 
the larger the misalignment angle 9t is, the larger is the 
magnitude of this contribution and, consequently, the 
larger is the increase in the wind velocity. 

Figure [6] shows line-radial cuts in the xz-plane of ra- 
dial velocity for the same instant shown in Fig. [S] Fig. |B^ 
presents radial cuts along the magnetic pole, i.e., along 
co-latitude 9 = 10° for T02 (where 9t = 10°), along co- 
latitude 9 = 20° for T03 (where 9t = 20°) and so on; 
Fig. IBb presents radial cuts along the magnetic equator, 
i.e., along co-latitude 9 = 100° for T02 (where 9t = 10°), 
along co-latitude 6* = 110° for T03 (where 9t = 20°) and 
so on. We note that along the magnetic pole, the curves 
of radial velocity for the tilted cases oscillate around 
the curve assigned for the aligned case (black solid line, 
9t — 0°). However, along the magnetic equator, the 
radial velocity increases as 9t gets larger, as explained 
in the previous paragraph. The wind radial velocity at 
- 75 ro is 227 km s'^ for case TOl, 234 km s'^ for T02, 
251 km s-i for T03, and 275 km s-^ for T04. Figures[6t 
and [6ji presents the same cuts along the rotational poles 
and equator, respectively. 

Figure [7| presents the inner portion of our simulation 
boxes at t = 240 h for cases TOl to T04. The magnetic 
field lines are represented by black lines, and the white 
line represents the contour-line where Br = 0. By fol- 
lowing the white line in Fig. [3 we note that the closed 
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Fig. 3. — Temporal evolution of the fiducial case T04 during one full stellar rotational period. The panels show 9 different instants, 
starting from t = ti = 216 h, and increasing in multiples of 1/8 Pq = 3 h, until t = ti + Pq = 240 h. We plot meridional cuts of the total 
wind velocity utot ■ The magnetic field lines are represented by black lines, and the white line is the contour-line of Br = 0. The x-axis 
ranges from to 10 ro and the z-axis ranges from —5 ro to 5 ro. 



magnetic field lines, as well as the open field lines, are 
not rigid, presenting a warped zone around the rotational 
equatorial plane of the star {z — 0). The amplitude of the 
oscillations gets larger as Of increases. The color maps 
show azimuthal velocity. There is not a significant vari- 
ation in the magnitude of u^p between the simulations, 
but the spatial profile of u^ is highly dependent on the 
configuration of the magnetic field, and thus on 9t. As 
we can see, the highest values of u^ are achieved inside 
the closed magnetic field lines (close to the star) and 
the rotating wind is forced to follow the same oscillation 
pattern of the magnetic field lines. In case of perfect 
alignment, the maximum azimuthal velocity happens in 
the equatorial plane of the star (Fig. [7^) . 
Two further simulations for different misalignment an- 



gles were performed: T05 (6't = 60°) and T06 {Ot = 90°). 
These simulations represent more extreme case of mis- 
alignment, being T06 the case where the axis of the mag- 
netic moment at the base of the coronal wind is perpen- 
dicular to the rotational axis. Both of them presents 
similar characteristics as the cases presented so far for 
9t < 30°, with an enhanced wind velocity though. For 
case TOG, the region of lower radial velocity remains in 
the region of closed magnetic field lines, but this region is 
now around the rotational poles (at co-latitudes 6* = 0°, 
180°). Figure [5] presents the 3D view of selected mag- 
netic field lines for cases T05 (Fig.|H]a) and T06 (Fig.|Hb), 
where it illustrates the inherent three-dimensional nature 
of our simulations. 
We selected case T04 to describe the wind character- 
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Fig. 4. — Three-dimensional view of the isosurface defined by 
Br = for case T04 after 8 rotations of the star (t = 192 h). The 
star is shown in red. 

istics. The periodic movement of the steUar magneto- 
sphere affects the entire wind structure, as can be seen 
in Fig. |9l where we present meridional cuts of the fol- 



lowing wind variables: Wtot) 
Meridional cuts of Ur and 



\Br\, \Be\, \B^ 



and Jtt 



can be found in Figs. [S] 



and [71 respectively. |i?r|, |i?e| and \B^\ present dipolar 
configuration at the base of the coronal wind, but for 
other radii, their solution is dependent on the resultant 
interaction of the magnetic field with the wind. Because 
of this interaction, the stellar magnetosphere acquires an 
azimuthal component for the magnetic field, which can 
be seen in Fig. I^Ji, presenting maximum intensity in the 
interface between closed and open field lines. The den- 
sity of the wind is not spherically symmetric, presenting 
higher densities around the Br = surface (Fig.|9j3). The 
total current density Jtot oc |V x B| is shown in Fig. |9f. 

3.3. The effect of a different 7 on the wind 

We now compare the results from simulations T04 and 
TO 7, where different values of 7 were adopted. The value 
of 7 influences the thermal acceleration of the wind. It 
defines the input of thermal energy and 7 obeys the re- 
lation p (X p'^ . Low values of 7 imply a proportionally 
large input of thermal energy in the wind, and conse- 
quently, high wind terminal velocities. It was shown in 
iVidotto et al.l |2009b) that when the wind is magnetized, 
the value chosen for 7 can alter the ratio between thermal 
and magnetic forces, and thus, accentuate the latitudinal 
dependence of the wind. 

In this work, we do not invoke the physical processes 
that may cause a larger input of energy in our models. 
Nevertheless, we study the effects of a smaller 7 in the 
wind of a magnetized star with oblique magnetic geome- 
tries. Figure [TU] presents the total velocity of the wind 
for both cases. By comparing simulations T04 (7 — 1.2) 
and T07 (7 — I.l), we find that the terminal velocity 
achieved by simulation T07 is around 22% larger than 
the one achieved in T04. The wind temperature profile 
for both cases is different, as expected: at the equato- 
rial plane, near 10 ro, for instance, T ~ 3.4 x 10^ K for 
case T04 and T ~ 5.5 x 10^ K for case T07. Further- 
more, the magnetic field configuration differs from both 



cases, with case T07 presenting a more compact zone of 
close d field lines. These re sults agree with our previous 
ones ([Vidotto et al.|[20"09ai rbl). which show that the heat- 
ing parameter is a important in defining the acceleration 
of the wind and the magnetic configuration around the 
star. 

3.4. The effect of different stellar rotational periods 

Here we explore the rotational effects on the dynam- 
ics of the wind, keeping dt — 30°. We select cases T04, 
T08, and T09 to perform this comparison, where the 
stellar rotational periods are 1, 3, and 0.5 d, respec- 
tively. These periods of rotation are in the lower range of 
observed periods for ^ yTTSs (e.g., iBouvier et al.l 119931: 
Choi fc HerbstI [ToOfil: iHerbst et al.l 120021 : IRebull et al.l 



20061) . as to explore the maximum effects on the wind 



Longer stellar rotational periods imply a wind that is 
less disturbed by the precession of the magnetic field. 
Figure [11] shows the meridional cuts of total wind ve- 
locities for Po = 0.5 d (Fig. [nji), Pq = 1 d (Fig. [Hb), 
and Pq = 3 d (Fig. ITTb ) at the same stellar rotational 
phase. The wind is more accelerated for lower Pq (i.e., 
larger rotational velocities), as a r esult of the coupling of 
magnetic fields and rotation fe.g-. lWeber fc Davislll967l : 
iBelcher fc MacGregod 119761 ). At the equatorial plane, 
wtot = 360 km s~^ for case T09, Utot = 250 km s~^ for 
case T04, and Utot — 180 km s^^ for case T08. 

3.5. The effect of a different I3q on the wind 

In IVidotto et al.l ()2009a Jbl , we have shown that the 
ratio between the thermal and magnetic energy densities 
at the base of the wind (/3o) is a decisive factor in defining 
the magnetic configuration of the wind, as well as its 
velocity distribution. To study how /3o infiuences the 
wind profile in the case of an oblique magnetic geometry, 
we compare simulations T04 and TIO, which have the 
same model parameters, except for the density at the 
base of the wind (and thus, different /3q). Figure[T2]shows 
meridional cuts of the total velocity of the wind plotted 
for both cases, as well as the magnetic field lines (black 
lines), and the surface Br = (white line). Both panels 
show a snapshot at the same rotational phas e of the star. 
As in the aligned case (jVidotto et al.l[20"09allbh . the wind 
is more accelerated for low /3o, where the magnetic energy 
density at the base of the wind is more important than 
the thermal energy density. The ratio of open to closed 
magnetic field lines are larger for case TIO, with lower 
/3o, showing that a faster wind is able to open the field 
lines more efficiently. 

3.6. The "wavelength" of the magneto spheric oscillation 

The typical length A of the oscillation of the isocon- 
tour oi Br = Q (white lines in Fig; s . 3, 7 , 9 - 12) can be 
estimated as in iLazarian fc Opherl ()2009( ) 

A ~ UcharPo = 0.062 ro liihT/''^^'" ' (14) 

where Uchar is the characteristic velocity of the plasma. 
For Po = 1 d, and different Ot (cases T02 to T04), we 
showed that the velocity of the wind increases with Ot- 
Therefore, from Eq. p^ . it is immediate to conclude that 
A will be larger for larger 9t. In fact, our estimates for 
the inner portion of the grid show that A ~ 13, 14, and 
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Fig. 5. — Meridional cuts of radial velocities u,- for (a) the aligned case with 8t = 0° (TOl) and cases with three different misalignment 
angles (b) dt = 10° (T02), (c) 8t = 20° (T03), and (d) 0t = 30° (T04). The panels are snapshots of the simulations after 10 stellar rotations 
(t = 240 h). In this figure, we present the meridional cut of the entire simulation box. 



16 ro for 9t — 10°, 20°, and 30°, respectively, which are 
consistent to the vahies measured from the simulations, 
Asim ~ 14, 16, and 17 tq. 

In the cases where 7 was compared (cases T04 and 
T07), we note that A is relatively larger for case TOT 
(7 = 1.1) than for case T04 (7 = 1.2). This is again due 
to the larger velocities of the wind for case T07. The 
value of A calculated by Eq. (fT4|) was ~ 16 and 20 ro for 
T04 and T07, respectively, while the measured values of 
Asim from the inner region of the simulations were ~ 17 
and 19 vq. 

In the cases where Pq was varied from 0.5 to 3 days 
(cases T04, T08, and T09), the velocity of the wind de- 
creased with the increase of Pa. Both facts have a differ- 
ent effect on Eq. ([T4| . However, the increase in Pq has 
proved to be more important in the increase of A than 
the effect provided by the increase in the wind velocity. 
We find A ~ 11, 16, and 37 tq for Pq = 0.5 d (T09), 
Po = 1 d (T04), and Pq = 3 d (T08), respectively, while 
the measured values were Asim — 12, 17, and 36 tq. 

Varying the ratio between the thermal and magnetic 
energy densities at the base of the wind (cases T04 



and TIO) has an important effect on the acceleration 
of the wind. Because case TIO presents a lower /3o than 
case T04, TIO has a larger characteristic velocity of the 
plasma. Ultimately, this increases A from ~ 16 ro (T04) 
to 30 To (TIO). The measured values of Asim from the 
inner region of the simulations were 17 and 31 ro. 

4. DISCUSSION OF STELLAR WIND RESULTS 

A MHD wind model provides solutions for the den- 
sity, velocity, and temperature profiles, along with the 
magnetic field configuration of the wind. Because the 
solution of the MHD equations is complex, some mod- 
els adopt several appro ximations. This is th e case of 
the Weber-Davis model () Weber fc Davis! [l967l ). first de- 
veloped for the solar wind and lat er on adopted to de- 
scribe winds of oth er stars (e.g., iPreusse et al.l 120051 : 
iLovelace et alJ 120081 among many other applications). 
The simplifications involved in the Weber-Davis model 
are: the model is axisymmetric and stationary; it con- 
siders an open, radial magnetic field that acquires an 
azimuthal component due to the rotation of the star; 
the wind solution is valid for the equatorial plane; it 
neglects meridional components of the magnetic and 
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velocity fields. Because only the open magnetic field 
lines contribute to angular momentum loss, the Weber- 
Davis model is expected to overestimate angular mo- 
mentum loss through a magnetized stellar wind, pre- 
senting shorter time-scales for stellar rotation brake than 
those model s that consider the exi stence of closed field 
line regions (jMestel &: Spruit! Il987l ). Because of its uni- 
dimensional characteristic, the solution of the Weber- 
Davis model can be easily integr ated. A deta i led de - 
scription of its solution is given in iPreusse et al.l ()2005D . 
Opposed to the Weber-Davis model, our model 
presents a multi-component corona, with the co-existence 
of open and closed field line regions. A latitudinally 
dependent velocity is observed, where the wind along 
the magnetic poles has a larger velocity than the wind 
along the equatorial regions. Details of the characteris- 
tics of the solutions of our wind model, in the context of 
aligned rotat ional axis and magnet ic dipolc moment, are 
presented in IVidotto et al.l (|2009bl ). Unfortunately, the 
three-dimensional nature of our model does not allow for 
analytical expressions of the solution of the variables of 
the wind. 



Different wind scenarios on the fram ework of the 
Weber-Dav is model were explo r ed b y IPreusse et al.l 
t2005i) and iHolzwarth fc Jardini (pOOl . The first one 
investigated the characteristics of the stellar wind for 
a sample of stars hosting close-in giant planets, while 
the second one used empirical data to constrain theo- 
retical wind s c enario s. De spite the wind solutions from 
IPreusse et all (|2005[ ) and IHolzwarth fc Jardind (|2007D 
being focused mainly on cool main-sequence stars, we 
compare the overall trend of our model in respect to 
these works. Mass-loss rates on such mod els are similar 
to th e solar wind value (10"^'* M© yr~^ IPreusse et al.l 
120051) and do not exceed ab out 10 times the solar value 
( Holzwarth fc Jardind l20Q7t). Wind terminal velocities 
found by IPreusse et al.l (|2005[ ) range from 310 km s~^ 
for an isothermal wind temperature of Tq = 5 x 10^ K 
to 760 km s~i for Tq = 2 x 10^ K. Compared to our 
models, mass-loss rates are about 6 orders of magnitude 
smaller, while velocities achieved are about the same or- 
der of magnitude. The difference in mass-loss rates is 
a consequence of the larger coronal densities adopted in 
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our models. 

To our knowledge, there are no measurements of mass- 
loss rates and wind velocities for WTTSs to compare our 
results with. The detections of mass-loss rates are based 
on the early stag e as a CTTS, when an accre tion disk is 
still present (e.g.jKuhilll964l: Edwards et al.ll2 003. 2006: 
IKwan et all I2007t IGomez de Castro fc Verdugg i2007i ). 
Based on these detections, mass-loss rates are of the or- 
der of lO"^'^ to 10~^ M0 yr~^ and wind terminal ve- 
locities ~ 400 km s^^. The parameters in our models 
were chosen to have values compatible to these ones. 
I Vidotto et al.l ()2009aD considered models with the den- 
sity at the base of the coronal wind spanning by two 
orders of magnitude, which resulted in mass-loss rates 
ranging between ~ 10~^ and 8 x 10~^ Mq yr~^. In the 
present paper, the main goal of iJ3] was to analyze the 
effects of the tilt angle on the wind. We therefore did 
not explore several values of base density and our mod- 
els present mass-loss rates of about 9 x 10~^ Mq yr~^. 
From Figures [3 [101 El and [12] we note that the wind 
terminal velocities obtained are ~ 350 to 500 km s^^. 



Kwan et al.l (pOOTJ ) 
WTTSs are simply 



suggest that, 
a scaled-up ve 



if 



the 
I of 



winds of 
the solar 



wind, WTTS winds should then be stronger than those 
of CTTSs, because X-ray emission from WTTSs are 
stronger than the emission from CTTSs. However, the 
winds of CTTSs are believed to be powered by accretion, 
which could be the reason why the wind traced by He l 
A10830 is not detected in WTTSs ()Kwan et al.|[2007[) . 
We would expect that when accretion ceases, the wind 
should become less strong. 

There is clearly a need of more observational constrains 
on the winds of WTTSs. In possession of that, we would 
be able to better constrain the parameters of our models. 

5. ON THE PLANET- WIND RECONNECTION 

Finding planets around young stars is currently ongo- 
ing, with two recent detections of mass ive giant plan- 
ets: one around a 5 Myr-old star (Lafrcni ere et al.ll201Clf) 
and one around a 12 Myr-old star (Lagra nge et al.ll2010D . 
The stellar wind is expected to directly influence the 
planet and its atmosphere, e.g., by changing the config- 
uration of the planet's magnetosphere, producing non- 
thermal planetary magnetospheric radio emissions, etc. 
So far, the few theoretical works investigating the influ- 
ence of the stellar wind on the magnetosphere of planets 
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respectively. 



Fig. 8. — 3D view of selected magnetic field lines for (a) 0t = 60° 
(T05) and (b) Ot = 90° (T06). 

were based on simplified treatments of the stellar winds, 
e.g., u si ng the Parker w ind model (jGrieBmeier et all 
l2007bllal:lLazio e t al."201(y). the Weber-Davis wind model 
(jPreusse et al.l 2005. 2006, 2007 ), assuming a s o lar- type 
stellar wind [iv ct al. 2004; iFarrell et al.l 120041 : iStevensI 
12005: Zarka 2007), or based on scali ngs for the mass- 
loss rates a nd wind terminal velocities (jGriefimeier et all 
12004 I2005D . The consideration of a realistic wind is cru- 
cial to determine how the interaction between the stellar 
wind and the magnetosphere of an extrasolar planet oc- 
cur. 

Our 3D, time-dependent MHD simulations of stellar 
winds of WTTSs provide a powerful tool to investigate 
the planet-wind interaction, as it allows us to consider 
the effects of a more realistic wind and obtain key in- 
sights on the detectability of radio emission from extra- 
solar planets. In this section, we estimate the reconnec- 
tion rate and power released when reconnection between 
a close-in magnetized giant planet and the stellar wind 
takes place. This estimate is performed for four different 
wind simulations: TOl, T02, T03, and T04, where the 
misalignment angle of the stellar rotation axis and its 
magnetic moment vector is 6t = 0°, 10°, 20°, and 30°, 



5.1. The Reconnection Rate 



In this section, we estimate the rate of reconnection 
between the stellar and planetary magnetic field lines. 
This is necessary to evaluate the planetary radio emission 

(gSXll). 

The magnetic field of the stellar wind has three com- 
ponents Bx, By, and Bz- B^ and By are parallel to the 
stellar rotational equatorial plane. Bz is perpendicular to 
the this plane. Considering a planet whose orbital plane 
coincides with the rotational equatorial plane of the star 
and considering that the planet's magnetic dipole mo- 
ment is aligned in the — z-direction, then when the mag- 
netic field of the stellar wind and the magnetic field of 
the planet are oriented anti-parallel to each other, mag- 
netic field line reconnection can occur. This results in 
transfer of energy from the stellar wind to the planet's 
magnetosphere. 

The reconnection rate is the amount of magnetic flux 
that reconnects per unit time per unit length of the re- 
connection line or, equivalently, the reconnection rate can 
be defined as the strength of the elect ric field parallel 
to th e reconnection merging line (e.g., iPriest fc Forbes! 
I2000r) . The reconnection line refers to the line where 
magnetic field lines reconnect. The rate at which re- 
connection between two different plasmas happens de- 
pends, among other things, on the ve locity of the inci- 
dent plasma on the reconnection site ()Parkeii[l973l) . In 
the idealized case when the magnetic fields of two iden- 
tical plasmas are exactly anti-parallel, the reconnection 
rate (or the generated electric field at the reconnec tion 
site) can be estimated as fe.g.. Borovskv et al.ll2008D 



E '. 



Vi-nB 



(15) 



where «;„ — Cva is the inflow velocity of the plasma 
in the reconnection site, va and B are the Alven speed 
and magnetic field of the ambient plasma, respectively, 
and c is the speed of light. The factor C — 1/ L is the 
dissipation region aspect ratio, i.e., a property of the 
geometry of the reconnection region which has a char- 
acteristic width / and a characteristic length L. When 
reconnection occurs between two plasmas with different 
characteristics, Eq. (|15p becomes more complex, taking 
into account the different magneti c field intensities and 
Alven speeds of the two plasmas (jCassak fc Shavl [20071 : 
IBorovskv erani2008l) 



E~C 



R3 R3 



1/2 



c \At:{Bz,2Pi+ Bz,iP2)(.Bz,i+ Bz^2 



(16) 



where the index "1" and "2" are used to distinguish both 
plasmas on the site of the interaction, B^ i and Bz^2 are 
oriented anti-parallel to each other, and p is the mass 
density. 

5.2. Planetary Radio Emission 

The solar wind interaction with the magnetic plan- 
ets of the Solar System (Earth, Jupiter, Saturn, Uranus, 
and Neptune) accelerates electrons that propagate along 
the planets magnetic field lines, producing electron cy- 
clotron radiation at radio wavelengths (|Zarka,1998, ). By 
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Fig. 9. — Meridional cuts of several wind variables for case T04 (6t = 30°). The panels are snapshots of the simulation at t = 240 h. 
Color-maps are for the (a) total wind velocity Mtot; (b) radial component of the magnetic field |-Br|; (c) 6-component of the magnetic field 
|i3g|; (d) azimuthal component of the magnetic field |iJ(^|; (e) mass density of the wind p; (f) total current density Jtot. All color-maps, 
except panel (a), are in logarithmic scale. 
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Fig. 10. — Meridional cuts of the total wind velocity «tot for 
cases with (a) 7 = 1.2 (T04) and (b) 7 = 1.1 (T07). Black and 
white lines have the same meaning as before. 

analogy to the magnetic planets in the Solar System, 
predictions have been made that magnetized extra-solar 
planets should produce c yclotron maser emission (e.g. , 
Farrell et all [1991 [200l IZarka et al.l [20011: jLazio et al 






20041 : ICxriefimeieret al.l [2004 120051 20n7ailbl : ISte^enf 
200a lZarkall2007t iJardine fc Cameronl l2008[ ). 

The planetary radio emission depends on the planet's 
magnetic field intensity^ and on the stellar wind power: 
in general, it implies that the stronger the stellar wind 
is, the more radio-luminous should the planet be. So 
far, such radio signatures from stars hosting hot Jupiters 
have not yet been detected, and one possible reason for 
that may be due to the lack of instrumental sensitivity 
in the appropriate f requency range of the observations 
(jBastian et al.|[2000D . The theoretical estimates on the 
radio flux emitted by extrasolar planets carry along a big 
uncertainty due to the fact that the stellar wind proper- 
ties are poorly constrained: iFarrell et al.l ()1999( ) showed 
that a variation by a factor of 2 in the wind velocity 
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Fig. 11. — Meridional cuts of total wind velocities utot at the 
same stellar rotational phase for (a) Pq = 0.5 d (T09), (b) Pq = 1 d 
(T04), and (c) Pq = 3 d (T08). Black and white lines have the same 
meaning as before. 



^ For pre dictions of pl a netary r adio emissioii fo r non-magnetized 
planet, see IZarka et all 11200 ij) : IZarkal I I2007I) : IGriefimeier et al.l 
ll2007al'l . 
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may change the level of radio power emission by a factor 
of 100, with more recent works suggesting that the ra- 
dio power ernission is proportional to the inciden t wind 
power (e.g., Stevens 2005; Zarka 2007; Griefimeie r et al.l 
l2007al. Therefore, the potential observation of radio 
emission from extrasolar planets strongly depends on the 
nature of the stellar wind. 

Based on our simulated stellar winds (cases TOl to 
T04), we estimate the planet's radio power. The elec- 
tric field generated in the interaction is calculated from 
Eq. (|16l) . where plasma 1 refers to the characteristics 
of the planet's magnetosphere, while plasma 2 refers to 
the local characteristics of the impacting stellar wind. 
Initially, we do not consider pile-up of the stellar wind 
magnetic field in the magnetosheath of the planet, but 
will do so in iJ5.2.3l 

To derive the characteristics of the planet's magneto- 
spheric plasma (plasma 1), we consider a hot Jupiter 
with a dipolar magnetic field aligned in the — z-direction, 
and magnetic intensity at the equator Bp = 50 G. The 
density of the planetary plasma pi is taken to be neg- 



ligible such that p2Sz,i > PiBz,2 in Eq. ^0^-^. We 
assume that the planet has the same radius as Jupiter 
Rp = Rjup ^ 0.05 tq. 

As we do not include the planet in our simulation, next, 
we analytically calculate the area of the planet that will 
interact with the stellar wind. 

5.2.1. The Size of the Planet's Magnetosphere 

The interaction of the planet's magnetosphere with 
the wind takes place at a distance tm from the cen- 
ter of the planet, where there is balance between the 
wind t otal pressure and the magnetic pressure of the 
pknet ("Lovelace et all l2008t iJardine fc CameronI 120081: 
rVidotto c t al. 2009c|'^ 



P2ul 



47r 



An 



(17) 



where Bz,i is the z-component of the magnetic field of 
the planet at the equatorial plane 



S.,i = 



BpRp 



' M 



(18) 



i?ll^2 is the parallel component of the stellar magnetic 
field to the boundary layer and U2 = (u^ — uk) is the 
relative velocity between the wind azimuthal velocity u,^ 
and the circular Keplerian velocity of the planet uk — 
(GAf*/r)i/2. Substituting Eq. ([HJ) in ^, we have 



Rp 



Bll2n 



{P2ui + By2ii) 



1/6 



(19) 



The size of the planet's magnetosphere tm depends on 
the local characteristics of the stellar wind, on the or- 
bital radius of the planet (through uk), and on the plan- 
etary magnetic field. Figure [T3b presents the size of 
the planet's magnetosphere if the stellar wind is given 
by cases TOl (aligned case, Ot = 0°), T02 [Ot = 10°), 
T03 [Ot = 20°), and T04 {Ot = 30°), calculated at 
t = 240 h, for a range of planetary orbital radius up 
to 12 To ^ 0.11 AU. We note that the planet's magne- 
tosphere becomes large as the misalignment angle 9t is 
smaller (i.e., the highest values of vm are found for the 
aligned case). This is a result of a lower wind total pres- 
sure ^(^2^2 + Bu 2/27r) as 9t gets smaller. For the range 
of orbital radii analyzed here, both kinetic and magnetic 
terms of the wind total pressure make significant contri- 
butions, except for small orbital radii (?" ^ 3 ro), where 
the magnetic term dominates. We see that for r < 2 rg, 
the magnetosphere of the planet has vanished, due to a 
large wind total pressure. As for the misaligned cases 
the wind impacting on the planet changes characteristics 
according to the stellar phase (see for instance Fig. [3]) , it 
is expected that the radius of the planet's magnetosphere 
will vary if the stellar rotational period is different from 
the orbital period of the planet. Considering a planet 
located at an orbital radius of r = 5 tq ~ 0.05 AU in 

^ In our simulations, for r < 12 ro, the ratio 1-62,2/^2,1 1 ranges 
from 0.2 to 9.3. For example, for case T02 at r ~ 5 ro ~'0.05 AU, 
-Bz,i ~ 6.1 G and -62,2 — —3.4 G. As -62,1 and -82,2 have approx- 
imately the same order of magnitude, this condition implies that 
Pi < P2 ~ 10~l^ g cm~3. 
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Fig. 13. — The magnetospheric radius tm [Eq. I|19| l] of a planet 
orbiting a star as given by simulations TOl, T02, T03, and T04. 
(a) Considering a given phase of the stellar rotation (4 = 240 h), 
and planet located in the x-axis. (b) Considering an orbital radius 
of r = 5 ro ~ 0.05 AU as function of stellar rotational phase 
(lines presented are merely symbol connectors and do not represent 
actual values of r^i). 



a circular orbit, for case T02, for instance, the variation 
in Tm from its maximum possible value to its minimum 
possible value is around 5%, while for cases T03 and T04, 
this variation is around 8% and 11%. Figure [TBb shows 
Tm as a function of the stellar phase of rotation for a 
planet located at r = 5 ro — 0.05 AU. 

Knowing the value of tm^ we can then calculate from 
Eq. (fT8| the value of the magnetospheric magnetic field 
of the planet B^^i that will interact with the stellar wind. 
Because the size of the planet's magnetosphere can in- 
crease/decrease depending on the incident wind, Bz,\ will 
present variations along the planetary orbit. The values 
of Bz 1 can range between 6.1 and 7.0 G for case T02, 6.8 
and 8.9 G for case T03, and 7.4 and 10.4 G for case T04, 
for a planet at an orbital radius r = 5 ro — 0.05 AU. 
Thus, for instance, for dt = 30°, along the planetary or- 
bit, the interacting planetary magnetic field varies up to a 
factor of 1.4. The more misaligned is the stellar rotation 



axis in relation to the stellar magnetic moment vector, 
more variation is expected in the magnetospheric radius 
of the planet and on the interacting planetary magnetic 
fields^,!. 

5.2.2. Estimate of the Planetary Radio Emission 

The power released with the reconnection event Prcc 
can be decomposed into a power released from the dis- 
sipation of kinetic energy carried by the stellar wind Pk 
and a power r eleased from the magnetic energy o f the 
wind Pb (;e.g.. lFarren"eran[T99l[Zarka et al.ll200l 

Prec = aPk + hPg , (20) 

where a and b are efficiency ratios and we assumed that 
Prec depends linearly on Pk and Pb- Observations of the 
Solar System can be either explained by {a = 1 x 10~^, 
h = 0) or \a = 0, 6 = 2 x IQ-^} ([Zarka 2007). In fact, 
iZarkal (J2007I ) argues that it is not possible to decide which 
incident power actually drives the radio power observed 
from the magnetic planets of the Solar System. If both 
incident powers contribute to the radio emission, this 
implies that the coefficients a and h need to satisfy the 
relation a/(l x 10"^) + 5/(2 x 10"^) = 1, as to match 
the observed radio power. However, we do not know 
if a and h should remain the same in other planetary 
systems. As we lack a better guess, we will adopt a and 
b as in the Solar System, i.e., either P^cc = 1 x lO^^P^ 
or Prcc = 2 X IQ-^Pb- 

The magnetic power Pb can be estimated as the Poynt- 
ing flux of the stellar wind imp acting on th e planetary 
magnetospheric cross-section 5* ()Zarkall2d07l ) 



Pr = 



ExB 

c — : • dS 

47r 



EB 



z,2 



47r 



' M ; 



(21) 



where the electric field E is given by Eq. (|T6)) . The con- 
stant C in Eq. (fT6|) is assumed to be C = l/L ~ 0.1, 
as derived by differe nt analytical and num erical methods 
(for a discussion, see lBorovskv et aLll2008[ ). This implies 
that we are assuming that the reconnection happens in 
a region with a characteristic width I = O.IL, where L is 
the radius of the planet (the characteristic length). The 
kinetic power Pk is 



Pk 



■dS 



P2ul-Kr'lj , 



(22) 



where pram — P2u\ is the wind kinetic ram pressure. 

According to our stellar wind models, throughout the 
region of investigation (r < 12 ro ci 0.11 AU), we find 
that Pk > Pb ( the same behavior is found in Figure 
11 of IZarkal l2007f l. except near the orbital radius where 
the planet co-rotates with the stellar wind and, thus, 
"2 = {Uip — uk) — and P^ ~ 0. Figure [Ml shows the 
estimated powers that are released for the wind cases 
TOl [Ot = 0°, solid lines) and T04 {Ot = 30°, dot-dashed 
lines) as a function of planetary orbital radius for a single 
rotation phase at t = 240 h. 

The consideration of several phases of rotation of the 
star is taken in Figure [151 which shows the estimated 
power that is released for the wind cases T02 {Ot = 10°, 
Fig. [Hi), T03 {Ot = 20°, Fig.[TniD), and T04 {Ot = 30°, 
Fig. [l5bV Figure [TS] presents the power released from 
the incident magnetic power of the wind (Prec = bPB)- 
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Orbilal radius (r„) 



Fig. 14. — Estimated dissipated power P^cc due to the interaction 
of a hot-Jupiter with the stellar wind for cases TOl (aligned, solid 
lines) and T04 (8t = 30°, dot-dashed lines) as a function of orbital 
radius of the planet, considered to lie in the x-axis. Red curves 
consider the power released from the incident kinetic power (Prcc = 
aPk), while black curves are from the incident magnetic power 
(^rec = bPg ) . The curves refer to the solutions for the stellar wind 
at t = 240 h. 

The shaded area hes between maximum and minim.um 
power that can be released in the interaction, depending 
on the characteristics of the incident wind (i.e., phase of 
rotation of the star). We see that the maximum emit- 
ted power gets progressively high as 9t increases, while 
the minimum released power gets progressively low as 
9t increases. This causes the shaded area to be larger 
for case T04 than for case T02 or T03. For a planet at 
orbital radius r = 5 ro ~ 0.05 AU, the ratio between 
maximum and minimum released power due to a varia- 
tion in the incident wind is a factor of 1.3 for case T02 
{Ot = 10°), 2.2 for case T03 [Ot = 20°), and 3.7 for case 
T04 {dt — 30°). Figure [15] also shows the released power 
for the aligned case (red solid line), showing that an in- 
clination between the axis of rotation of the star and the 
surface magnetic moment vector can contribute for an 
increase in the emitted power, depending on the phase 
of rotation of the star. 

Part of this released energy can be used to accelerate 
electr ons, generating radio emission (jJardine fc CameronI 
[2?)n8h 



^radio "~ ^^rc 



(23) 



The efficiency rj in the conversion of Prec hito radio 
emission Pradio depends on the details of the physi- 
cal processes that generate the radio emission (e.g., on 
the cyclotron-maser instability). Assuming t] = 10% 
([jardine & Cameron 2008), the radio power emitted 
from a Jupiter-like planet orbiting at a distance r = 
5 ro ~ 0.05 AU is P,.adio -- 10^^ W. The radio power 
is also shown in Fig. [15] (vertical axes on the right). 
A time-depe ndent radio em ission has also been esti- 
mated by Fa res et all (|2010f l. For the magnetic plan- 
ets of the Solar System Pradio ~ lO^'^ W (Neptune) 
to ^ lO"'^^'^ W (Jupiter), which means that, for our as- 
sumed giant planet orbiting our fictitious star, the radio 
power released is almost ^ 5 orders of magnitude larger 
than for Jupiter. This result suggests that stellar winds 
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Fig. 15. — Estimated dissipated power Prcc = bPs due to the 
interaction of a hot-Jupiter with the stellar wind for cases (a) T02 
(et = 10°), (b) T03 {et = 20°), and (c) T04 {9t = 30°). The 
shaded area lies between maximum and minimum power that can 
be released in the interaction, according to the rotational phase of 
the star. The red curve is for the aligned case TOl {0t = 0° ). T he 
emitted radio power, assuming a conversion given by Eq. 1 )2311 . is 
shown in the vertical axes on the right of each plot. 
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from pre-main sequence young stars have the potential 
to generate stronger planetary radio emission than the 
solar wind. Our results are in accordance to previous 
works developed on the framework of st ellar winds of 
stars at t he early ma in-sequence phase (jStevensI 120051 : 



stars at t he early ma ir 
iGricfimcie ret al.ll2005f ) . 



Table [2] presents a summary of the results obtained for 
cases TOl to T04. The properties of the wind and of the 
reconnection site are described at r — 5 tq. 

5.2.3. Pile-up of the Stellar Wind Magnetic Field 

Around the Earth, the magnetic field of the solar wind 
piles-up in the magnetosheath. This causes the magnetic 
field strength to enhance in the reconnection site, and 
ultimately, increases the electric field E. By analogy, 
here we consider what would happen if magnetic field 
pile- up is considered in the planet's magnetosheath. 

If we consider that the stellar wind is supersonic, a bow 
shock forms and the wind is deflected around the planet 
magnetosphere. For cases TOl to T04, the wind becomes 
supersonic at r ~ 4 rO at the rotational equatorial plane. 
Downstream the shock, the field strength i?||.2, the ve- 
locity U2, and the density p2 that appear on Eqs. (|16|) . 
(fT9|) . and ([21]) will be given by shock conditions instead 
of arising directly from our stellar wind model. In this 
case, we use the Rankine-Hugoniot jump conditions to 

determine the magnetic field intensity B^f^, density P2 , 



and velocity Mj on the magnetosheath ( Qpheii [20100 



Hence, except for Eqs. ([T6|) . p9|) . and (|2T|) . the equa- 
tions presented in this section remain the same. 
Assuming a perpendicular shock with strength S — 

P2 / P2 — 4 (i.e., the density in the magnetosheath is 
four times the value of the stellar wind density), the 
Rankine-Hugoniot jump conditions state that the mag- 
netic field at the magnetosheath is ^ 4 times higher 
than the local value of the wind magnetic field (i.e., 

4_B|| 2), while the velocity drops by a factor ^ 4 



B 



11,2 



(i.e., M^ 



U2/4). As a consequence, the size 



M 



of 



the planet's magnetosphere [Eq. p9)) ] diminishes with re- 
spect to rM when pile-up of the wind field lines was not 

considered. Within r < 12 ro, the ratio r^^ /'''m ^ 0.68 
for case TOl, and becomes slightly more significant for 

higher 9t, presenting r^ /vm ^ 0.63 for case T04. De- 

(s) 

spite the decrease in r\/, the increase in the reconnec- 
tion rate E''^'*^ [Eq. (|16p ] is such that the dissipated power 

^rcc — bPg [Eq. ((2T|) ] increases when compared to the 
case when pile-up of the wind field lines was not consid- 
ered. For case TOl, Prcc increases 4.8 - 6.8 times the 
values of Prcc presented in Fig. [151 depending on the or- 
bital radius of the planet. For case T02, this increase 
ranges between 4.5 - 7.2, 4.5 - 8.4 for case T03, and 4.5 
- 10.5 for case T04, where these ranges depend now on 
both the location of the planet, as well as on the rota- 
tional phase of the star. This shows that for a supersonic 
stellar wind, the consideration of a perpendicular shock 
can increase the dissipated power due to the interaction 
of a hot- Jupiter with the stellar wind. 

5.3. On the Detectability of Planetary Radio Emission 



The detection of planetary radio emission depends on 
several factors, such as, on the distance d to the extra- 
solar system, on whether the conical beam of the cy- 
clotron emission is di rected towards us, and on the emis- 
sion bandwidth A/ ([Farreh et al.lll999[ ). The radio flux 
that we detect on Earth is given by 



*&radio — 



-'radio 

d^wAf ■ 



(24) 



where w is the solid angle of the conical emission beam. 
The stellar wind ultimately controls the incident power 
on the planet, while the planet's characteristics control 
the frequency of the cyclotron emission /c, and thus 
the emission bandwidth assumed to be A/ = 0.5/c 
([Farrell et al.|[r999.) .'^ For our fictitious planet, the as- 
sumed magnetic field at the pole is 100 G (maximum 
field strength), which emits at (maximum) fc — 2.8B — 
280 MHz (b given in G and fc in MHz), with a band- 
width of A/ = 140 MHz. If our star is at a distance 
d ^ 10 pc , using our estimated Pradio ^ 10^^ W obtained 
in i)5.2.2i we find that the radio flux detected at Earth 
would be $radio ~ 7.5/w mJy. For a spherical emission, 
w = An, and the detected flux is <i>radio ~ 0.6 mJy, while 
for a hollow-cone beamed emission with a conical aper- 
ture of 45°, w ^ 1.8 sr, and ^radio ^ 4 mJy. 

One of the possibilities of the yet unsuccessful ra- 
dio detections may be due to a probably small plane- 
tary field, which produces cyclotron emission in a low- 
frequency range, where in strumental sensitivity is still 
poor (Ba stian et al.l l2000t ). Low- frequency detectors, 
such as LOFAR, might be able to detect emission of few 
mJy at a frequency range of 10 to 240 MHz in the fu- 
ture ([FarreU et al.|[200llGnefimeier et al.ll2007aD . Our 
hypothetical planet, for instance, could be observable by 
LOFAR. 

Throughout ^ we have assumed a planet with equato- 
rial magnetic field of 50 G, which is about 6 times larger 
compared to Jupiter's mag netic field of 8 G at the equa- 
tor (or ~ 16 G at the pole, iConnernev et al.lll998f) . If we 
assume a planetary magnetic field intensity as Jupiter's, 
i.e., Bp = 8 G at the equator, the power Pioc = bPs 
released from the reconnection between planetary and 
stellar wind field lines at r = 5 ro — 0.05 AU will be 
a factor of 2.7 to 3.4 smaller than the estimates per- 
formed in H5.2.2l for Bp = 50 G. Also, because the plane- 
tary magnetic field is smaller, the magnetosphere of the 
planet will face a strong wind pressure and will vanish 
for r < A ro (in this case, tm — Rp)- For a value of 
Prcc — 3 X 10^^ erg s^^, and adopting the same efficiency 
rj = 10% for the conversion of the released energy into 
radio power, the radio flux of the planet arriving at Earth 
(adopting d = 10 pc) would be Cradle = 14/w mJy, at a 
bandwidth of A/ = 22.4 MHz, for a given solid angle w of 
the conical emission beam. Comparing to Bp = 50 G, the 
drop on the radio power caused by a smaller planetary 
magnetic field (Bp = 8 G) is more than compensated by 
a smaller A/ leading to a radio fiux that is about twice 
the value obtained for Bp = 50 G. 



* Other authors (e.g., GricB meier et ani2007al ) adopt a higher 
value of A/ ~ 0.9 — 1.0/c. If this is the case, the emission band- 
width we obtain will have to be multiplied by a factor of 1.8 - 2.0, 
and the resultant radio flux divided by the same factor. 
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THE INFLUENCE OF THE WIND ON PLANET 
MIGRATION 



In IVidotto et al.l (j2009af) , we investigated the action 
of magnetic torques from the stellar wind acting on the 
planet and whether such torques were able to signifi- 
cantly remove orbital angular momentum of the planet, 
causing planetary migration. The idea is that the wind 
exerts a pressure ptot = \{p2u\ + B?,^/'^^^) on the area 



ir 



This force 



of the planet's magnetosphere Acs — irrj^ 
Ptot^off will produce a torqu e -^ at the planetary or- 
bital radius r ([Lovelace et al.l [20081 ) 



dLr, 



dt 



{Vtot Acs)r 



(25) 



where Lp = MpVKr is the orbital angular momentum of 
the planet, and Mp is the mass of the planet. A change 
in the planet's angular momentum also leads to 



dLo 



dt 



1,. dr 
2 '' ^ dt 



1,. r 

-MpVK — 



(26) 



where r^, is the ti me-scale for an a ppreciable radial mo- 
tion of the planet ([Lin et al.lll"996[) . From Eqs. ^ and 
([25[) . we can estimate such time-scale 



1 MpVK 

2ptotAcs 



(27) 



Within the range of parameters adopted in the simu- 
lations performed in IVidotto et al.l ([2009a|) , we showed 
that the stellar winds of WTTSs were not expected to 
have strong influence on the migration of close-in gi- 
ant planets. The winds analyzed in that case assumed 
that the axis of rotation of the star and the stellar mag- 
netic dipole mome nt were aligned . One a spect that was 
not investigated in IVidotto et al.l ([2009aD is the effect of 
a tilted magnetic moment with respect to the rotation 
axis. Using the results of the simulations presented in 
the present paper, we thus compare the time-scales r^, 
obtained for the aligned case (TOl) and the misaligned 
case T04 {9t = 30°). We consider a planet with the same 
mass and radius as Jupiter, and the magnetic field at 
the equator of Sp = 50 G. This comparison is shown in 
Figure [161 where we note that an inclination of the stel- 
lar magnetic field acts to reduce r^ when compared to 
the aligned case. The cases with intermediate tilt angles 
investigated (i.e., 6t — 10°, 20°) results in time-scales 
Tw{r) that lie between the solid line for case TOl and the 
dot-dashed line for case T04. Figure [TCI illustrates one 
single rotational phase of the star. For other phases of 
rotation of the star, we observe the same behavior: r^, 
calculated for the misaligned cases is smaller than for the 
aligned case. 

The last column of Table [5] presents the time-scale t^ 
calculated for cases TOl to T04 for a planet at r = 5 tq. 
Compared to the aligned case, where r^, ~ 100 Myr, case 
T04 {9t — 30°), for example, shows considerably smaller 
time-scales ranging from r^, ^ 40 to 70 Myr. We expect 
that larger misalignment angles 9t > 30° or other effects, 
such as an increase in the wind coronal b ase density or 
magne tic field intensity (as discussed in IVidotto et al.l 
l2009a[) . could reduce r^,. The time-scales derived here 
seem to be larger (and therefore less important) than 




Orbital radius (r ) 



Fig. 16. — Time-scale t^ for an appreciable radial motion of the 
planet. Solid line is for case TOl (aligned case) and dot-dashed line 
for case T04 {6t = 30°). The planet is assumed to lie in the x-axis. 
The curves refer to the solutions for the stellar wind at t = 240 h. 

those estimated by other processes, such as from the in- 
teractio n of the protoplanet with the disk wherein it was 
formed (IPapaloizou fc Terauem|[2006l ). However, as sug- 
gested bv lLanzal ([20101 ). the removal of planetary orbital 
angular momentum by the stellar wind may be important 
for synchronizing stellar rotation with the orbital motion 
of their planets during the pro-main sequence phase. 

7. CONCLUSION 

We have presented simulations of magnetized stellar 
winds where the surface stellar magnetic moment is tilted 
with respect to the axis of rotation of the star. Such 
configuration requires a fully 3D approach, as the system 
does not present axisymmetry. By adopting a dipolar 
surface distribution of magnetic flux, we showed that the 
interaction of magnetic field lines and the wind leads to 
a periodic final solution, with the same rotational period 
as the star. The final magnetic field configuration of the 
stellar magnetosphere presents an oscillatory pattern. 

By varying several parameters of the simulations, we 
explored the effects of the misalignment angle dt , stellar 
period of rotation Pq, heating index 7, and plasma- /3 
parameter at the magnetic pole of the star /3o in the 
final periodic solution of our simulations. We showed 
that the increase in 6t or the decrease in Pq lead to a 
more accelerated wind. The same is true if 7 or /3q are 
decreased, as already demon strated in the axisym metric 
cases of our previous paper ([Vidotto et al.|[20"09aD . 

We estimated the power released in the interaction be- 
tween a close-in giant planet and the stellar wind. If the 
planet and wind are magnetized, the interaction results 
in reconnection process, which releases energy that can 
be used to accelerate electrons. These electrons prop- 
agate along the planet's magnetic field, producing cy- 
clotron radiation at radio wavelengths. This calculation 
is motivated by radio emission observed in the magnetic 
planets of the Solar System (Earth, Neptune, Saturn, 
Uranus, and Jupiter). We showed that the intensity of 
the radio emission varies, as the wind impacting on the 
planet changes according to the stellar phase of rotation. 
If radio emission from a planet orbiting a star with mis- 
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aligned rotation axis and magnetic axis is to be detected, 
we showed here that the radio power will present a larger 
temporal variation for higher 0t . 

Our estimates show that the radio power emitted by 
the fictitious extra-solar planet orbiting our star at ~ 
0.05 AU can be 5 orders of magnitude larger than the 
non-thermal radio power emitted by Jupiter. This sug- 
gests that the stellar wind from a young star has the po- 
tential to generate strong planetary radio emission, which 
could be detected by LOFAR. 

As a final point, w e answered the question posed in 
iVidotto et al.l (|2009a|) about whether winds from mis- 



aligned stellar magnetospheres c ould cause a significan t 
effect on planetary migration. In IVidotto et al.l (|2009aD . 
only the case of winds from stars where the rotation axis 
and the surface dipolar magnetic moment were aligned 



was considered. Compared to the aligned case, we 
showed here that the time-scale for an appreciable radial 
motion of the planet is shorter for larger misalignment 
angles. 
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TABLE 2 
Summary of the local characteristics of the wind and planet at the reconnection site for a planet at orbital radius r = 5 ro ~ 0.05 AU. The 

COLUMNS represent, RESPECTIVELY: THE SIMULATION NAME, THE LOCAL WIND DENSITY p2, THE RELATIVE VELOCITY U2 = (u^ — Ujf), THE MASS-LOSS RATE PER 

UNIT SOLID ANGLE pUrT^ , THE LOCAL Z-COMPONENT OF THE WIND Bz,2 AND OF THE MAGNETOSPHERE OF THE PLANET Bz,l, THE E ELECTRIC FIELD (OR 

RECONNECTION RATE) AT THE RECONNECTION SITE, THE SIZE OF THE PLANETARY MAGNETOSPHERE rju. THE POWERS RELEASED IN THE INTERACTION ADOPTING 

Pj.^^ = aPfc AND Proc = bPg, AND TIME-SCALE FOR PLANET MIGRATION T„ (JBJ. RANGE OF NUMBERS ARE DUE TO DIFFERENT PHASES OF ROTATION OF THE STAR. 

Name p2 "2 pUrV^ -Bz.2 Bz.l E Tm Prec = aPk Prec = bPB Tw 

(10^13 g cm-3) (kms-i) (IO-^'Mq yr^^) (G) (G) (10"^ stV cm-i) (Pp) (lO^^ g;.g g-l) (lo^^ erg s^i) (Myr) 
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